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Objectives of the Cooperative Agreement 


There were two principal objectives of the cooperative agreement between NASA and the 
University of Denver. The first goal was to modify the design of the ER-2 condensation 
nuclei counter (CNC) so that the effective lower detection limit would be improved at 
high altitudes. This improvement was sought because, in the instrument used prior to 
1993, diffusion losses prevented the smallest detectable particles from reaching the 
detection volume of the instrument during operation at low pressure. Therefore, in spite 
of the sensor’s ability to detect particles as small as 0.008 microns in diameter, many of 
these particles were lost in transport to the sensing region and were not counted. 

Most of the particles emitted by aircraft are smaller than 0. 1 pm in diameter. At the start 
date of this work, May 1990, continuous sizing techniques available on the ER-2 were 
only capable of detecting particles larger than 0. 17 micron. Thus, the second objective of 
this work was to evaluate candidate sizing techniques in an effort to gain additional 
information concerning the size of particles emitted by aircraft. 

Accomplishments during this Cooperative Agreement 

Redesign of the ER-2 CNC . The design of the ER-2 CNC was evaluated under this grant 
and revisions were made in the design which would permit a larger fraction of small 
particles to reach the detection volume of the instrument during low pressure operation. 
These design changes were implemented in the construction of the ER-2 CNC II which 
was used in the 1993 Stratospheric Photochemistry, Aerosol and Dynamics Expedition 
and in ASHOE/MAESA flights which are currently under way. Figure 1 shows the design 
changes incorporated into the ER-2 CNC as a result of this work. The redesigned 
instrument, the ER-2 CNC II, was constructed with support from another grant. The 
response of the current, redesigned CNC is less dependent on pressure than that of the 
previous instrument. 

Evaluation of sizing techniques. A number of techniques were evaluated for determining 
the sizes of particles smaller than 0.1 micron in diameter. These techniques included a 
differential mobility analyzer (DMA) and a diffusion battery (DB). In addition, the 
possibility of extending the lower size range of the optical particle counter was also 
considered. The object was to provide accurate size distributions of particles emitted from 
aircraft flying in the stratosphere. 

Evaluation of the DMA technique 

Technical reports generated in the study of the DMA are attached in Appendices 1 and 2. 
A summary of the results is given here. 

The DMA technique promises high resolution measurements of aerosol size distributions. 
The technique consists of exposing the aerosol to a bipolar ion atmosphere which will 



neutralize most of the particles but leave a certain number with a single charge, fewer with 
a double charge, etc. After this conditioning, the aerosol is introduced into the DMA 
where particles are exposed to known flow and electric fields. These fields are crossed so 
that only particles within a narrow range of electrical mobilities can reach the exit region. 
Thus, ignoring the effects of multiple charging, particles within a very narrow size range- 
virtually monodisperse-exit the DMA. These particles are counted by a CNC. By 
varying the voltage and/or the flow, a spectrum particle concentration vs. particle mobility 
may be measured. Accounting for multiple charging effects, this spectrum can be 
converted to a size distribution for particles larger than the CNC detection limit (about 
0.008 pm diameter). Thus, this technique hold the promise of high-resolution size 
distribution measurements of particles in the size range likely to be produced by jet 
engines. 

The analysis of the DMA alternative in Appendices 1 and 2 shows that there are two 
fundamental issues which remain unresolved for a stratospheric DMA. The first issue 
involves the effects of Brownian diffusion on the transfer function of the DMA. 
Discrepancies between the predicted and actual transfer function for very small particles- 
those most affected by diffusion— exist in the literature. These discrepancies suggests that 
the effects of diffusion are not precisely predictable or experimentally controllable. Since 
diffusion is increasingly important as pressure decreases, the response of a stratospheric 
DMA may not be folly predictable. The magnitude of this effect has not been ascertained. 

A second problem associated with the operation of the DMA at stratospheric pressures 
involves the charging of aerosol at low pressure (Appendix 2). The recovery of a size 
distribution from a DMA depends critically on the distribution of charges on aerosol 
particles resulting from the ion preconditioning process. Although ion properties and 
aerosol chargeing have been studied at atmospheric pressure, there is little literature on 
these processes at low pressure. It may be possible to use the preexisting ambient charge 
distribution and avoid the the artificial preconditioning problem. At the current time, the 
knowledge of stratospheric ion concentrations, ion-attachment coefficients at low 
pressure, and the charge distribution on particles found in aircraft exhaust is insufficient to 
evaluate the effectiveness of a stratospheric DMA. 

A third problem involving a stratospheric DMA, one not discussed in the appendix, is 
associated with the absolute concentration of particles available in the stratosphere. Only 
a small fraction of the particles upstream of the DMA are charged in the preconditioning 
process. Since the concentration of stratospheric particles is quite small, obtaining a 
reasonably accurate size distribution measurement in the stratosphere would require 
extremely long sampling times to obtain a statistically significant signal. The aerosol being 
measured would have to remain relatively constant over the long measuring time. Even in 
the plume of an aircraft, the large spatial variations in particle concentrations would make 
complete particle size measurements with a DMA impractical. 

Given the uncertainties in diffusion losses within a DMA at low pressure, the unknown 
distribution of charges on stratospheric particles, and the low concentrations of these 



particles that imply very long sampling times, using a DMA for measuring stratospheric 
aerosol size distributions appears to be impractical. Some of these uncertainties may be 
reolved with future laboratory measurements and modelling studies; however, the problem 
of long sampling times will probably require fundamental changes in DMA design and may 
remain unsolved. 

Evaluation of diffusion battery technique 

Diffusion batteries have been used for decades to evaluate the size distribution of small 
particles in the atmosphere. In this class of instruments, particles are separated by passing 
the aerosol through a “battery” of screens or tube arrays. Particles with high diffusional 
mobilities (i.e., small particles) diffuse to the walls of the battery and are removed from the 
airstream. The remaining particles are counted by a CNC. 

There are several problems associated with determining particle size distributions from 
diffusion batterie measurements. The main difficulty lies in the nature of the response 
function and the difficulty in recovering accurate particle size distributions from the 
response function. These difficulties are discussed in more detail in Appendix 3. Briefly, 
the diffusion battery does not provide a sharp discrimination of particle size, rather a very 
gradual one. The measured parameters are strongly convoluted with the response 
function and must be deconvolved, or inverted, to obtain the original distribution. 

Because this convolution is of the form of a Fredholm integral equation of the first kind, 
there are an infinity of possible solutions (particle size distributions) that could produce 
the measured response. Various techniques are availabe that constrain the equation and 
produce the “best” solution, while another reports all physically possible solutions, giving 
a wide range of particle size distributions. When measurement errors are added to the 
inversion, the recovery of a particle size distributions becomes very difficult. Because of 
these problems, D^As afd diffusion batteries are rarely used now in cases where a DMA 
can be applied. 

Because of the fundamental problems associated with inverting diffusion battery data and 
recovering a size distribution, it is unlikely the diffusion batteries can be used to evaluate 
the details of the size distribution of small particles in the stratosphere. However, the 
technique shows some limited promise in providing rough estimates of some important 
parameters— such as mean particle size-in the diameter range from 0.008 to 0.1 pm. 

Redesign of the PMS aerosol spectrometer 

In response to requests made by Dr. Guy Ferry of NASA Ames Research Center, Particle 
Measuring Systems Inc. (PMS, Boulder, CO) offered to build an optical particle counter 
capable of sizing particles as small as 0.05 pm in diameter. Since the instrument heats 
particles before measuring them and thus reduces their size, a sulfuric acid particle having 
a diameter of 0.05 pm at the laser would have a diameter approaching 0.07 pm under 
ambient temperatures and water vapor mixing ratios. This instrument, the focused-cavity 
aerosol spectrometer, or FCAS, represents a considerable advance on the previous model, 



the active-scattering aerosol spectrometer. NASA has purchased the upgraded optical 
particle counter. The FCAS flew successfully on the AASE-II and SPADE missions 
before being replaced with a model with an improved laser and new pulse analysis 
hardware. The FCAS instruments have collected data that indicate the presence of a mode 
of small particles when the CNC counts high concentrations of particles, although this 
mode is only partially resolved by the spectrometer measurements. 


Summary 

The ER-2 CNC has been modified to reduce the diffusive losses of particles within the 
instrument. These changes have been successful in improving the counting efficiency of 
small particles at low pressures. 

Two techniques for measuring the size distributions of particles with diameters < 0.17 pm 
have been evaluated. Both of these methods, the DMA and the diffusion battery, have 
fundamental problems that limit their usefulness for stratospheric applications. We cannot 
recommend either for this application. Newly developed, alternative methods for 
measuring small particles include inertial separation with a low-loss critical orifice and 
thin-plate impactor device. This technique is now used to collect particles in the multi- 
sample aerosol collector housed in the ER-2 CNC-II, and shows some promise for particle 
size measurement when coupled with a CNC as a counting device. 

The modified FCAS can determine the size distribution of particles with ambient diameters 
as small as about 0.07 pm. Data from this instrument indicate the presence of a nuclei 
mode when the CNC-II indicates high concentrations of particles, but cannot resolve 
important parameters of the distribution. 

Bibliography. 

No publications resulted from these engineering and feasibility studies. 

Figure Captions: 

Figure 1 . Assembly A shows the design of the ER-2 CNC. In this instrument, the sample 
was drawn through the sample flow meter, Al, up the injector in the vertical saturator, 

A2, and into the condenser. In the revised design, Assembly B, the flow meter and 
injector are combined which leads to a shorter length from the inlet to the saturator. Thus 
fewer particles are lost by diffusion. 

Figure 2. Modelled response function of a tube-array diffusion battery under stratospheric 
conditions. 
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Appendix 1 


STRATOSPHERIC DMA 


NOMENCLATURE 

Dimensional variables: 
t = time 

t = radial coordinate 
z— axial coordinate 
ft = flow velocity at r,z 
E = electric field at r,z at time t 
V = analyzer voltage at time t 

Dimensional parameters: 

Z p = particle electric mobility 

D = particle diffusivity 

R l = radius of central rod 

R 7 = inner radius of outer cylinder 

D h = 2{Ri-R\) = hydraulic diameter 

L = axial distance between aerosol entrance and exit slits 

Az a = axial width of aerosol entrance slit 

Azs = axial width of aerosol exit slit 

Q & = aerosol entrance Cow 

Q s = aerosol exit flow 

Q c = sheath entrance flow 

Q m = sheath exit flow 

Qt = Qa.+ Qc = Qs+Qn = total flow 

U = 0 t /7r(f? 2 2 -#i 2 ) = mea n flow velocity 

Vq = analyzer voltage at t = 0 

r = analyzer voltage scan time constant 

Z p * = (C? m +<?c)*log(Ri/R2)/47rLU= centroid of steady-state transfer function at V 
Z p o* = Zp’{ V= Uq) 

Dimensionless variables: 

1 = t/r 



0 ;=(r/fl 2 ) 2 
z = zj Dh 

Dimensionless parameters: 

7 = (R]j R7) 2 
K = LR2/{R7 2 -Rl 7 ) 

Ka. = Az a /T 

Ks = AZs/L 

/? = ( <2s"b Qa)/{ Qin+ Qc) 

5= {Qs-Qz)KQs+Qz) 

G = geometric flow factor 
D= 4irLDl{Q m +Qc) 
a 2 = G- D 

a * 2 = a 2 (^ p =^p*) = G-D* 

Re = pUDh/n = Reynolds number' 
5= L/rU 
~Zp — Zpj Gpo* 


INTRODUCTION 

It is proposed to develop a new particle sizing system to be operated aboard the 
ER-2 in the stratosphere. One of the possible measurement systems consists of a 
differential mobility analyzer (DMA) operated upstream of a condensation nucleus counter 
(CNC). The purpose of the DMA is to extract aerosol within a known narrow size range 
and present it to the CNC for concentration measurement. The response of the DMA is 
well understood in its standard steady-state mode (Knutson and Whitby, 1975a). 

However, for this project it proposed to operate in the scanning mode described by Wang 
and Flagen (1990). Use of the scanning mode may present several problems depending on 
the rate of scan. Some of these problems may be confounded by particle Brownian 
diffusion which is significant for the proposed application. Aerosol charging is also a major 
area of uncertainty. 

In the discussion below, the DMA scanning mode is first analyzed neglecting 
diffusion. Diffusion effects are discussed in the following section. Aerosol charging is 
discussed in a separate paper. 
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SCANNING MODE AND CHARACTERI STIC TIMES 

Strictly laminar flow may be assumed within the DMA classification section itself. 

In general, laminar flow may be assumed from some indeterminate point upstream of the 
DMA all the way to the CNC viewing volume. There may be some mixing at certain 
points in the transport lines, such as just inside the DMA aerosol exit slit, which will affect 
the following analysis. But for the most part, particles may be assumed to follow flow 
streamlines in the absence of diffusion and the DMA electric field. 

In analyzing the transit of a particle through the proposed measurement system, 
there are several key locations along its trajectory from the sample point to the detector. 

The first is where it enters the sampling system at the tip of the external sample probe, 
designated as p. Well mixed flow is assumed from p to b, the point where laminar flow 
begins upstream of the DMA. The DMA aerosol entrance and exit slits are designated as a 
and s, respectively, the CNC detection volume as d. Let r m represent the average particle 
transit time through the well mixed zone from p to b. Let r a , r x and r s represent the 
average particle transit times from Ho a, a to 5, and 5 to d , respectively. Let £p, tb, k, k, 
and t d represent the times at which a particular particle passes the respective points. An 
inherent assumption for the well mixed zone is that all particle transit times are the same, 
i.e. f b -£ p =r ra . But for the laminar flow regions the transit times for a particular particle 
depend on which streamline it is on. Because of slow flow near walls, the variation of 
transit times through the laminar zones upstream and downstream of the DMA is on the 
order of the mean transit times, r a and r s . The variation of transit time through the DMA 
itself is significantly less than r x . Let r represent the period over which the voltage of the 
scanning DMA changes by a significant amount, and rf the period over which the ambient 
aerosol size distribution changes by a significant amount. The CNC sample interval is 
designated as r d . 

Even when operated in the steady-state mode, the size range extracted by the DMA 
has a finite width determined by the ratio of aerosol to sheath flows. When operated in the 
scanning mode there are several factors acting to broaden the particle size range presented 
to the CNC at any given instant. 

Particles arriving at the CNC at a given instant have exited the DMA at point s 
over a range of earlier times. If r s <<r does not hold, then the mean particle size exiting 
the DMA has changed significantly over this time range resulting in a broadening of the 
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size range presented to the CNC at that instant. 

If t x «t, then the DMA voltage may be treated as constant over the period of a 
single transit of the DMA classification section and there is no broadening of the extraction 
range for this part of the particle trajectory. However, as the scan rate is increased (r 
decreased) the DMA transfer function is increasingly warped from its basic triangular 
shape and becomes broadened. (See results of SCAN simulation program.) 

If T d <<r does not hold, then again the particle size exiting the DMA has changed 
significantly over the sample interval resulting in broadening of the effective average 
extraction range over the interval. 

For normal stratospheric sampling r f may be assumed to be large compared to the 
sampling system transit time. This condition may be strained somewhat for the High 
Speed Research Program when flying through a fresh jet plume. If the laminar flow zone 
transit time is not small compared to r f then particles arriving at the CNC at a given 
instant have been sampled from the ambient over a range of times during which the 
ambient size distribution has changed significant!.) . 

In general, if r and r f are large compared to r a , r x , r s and rj then the 
measurements can be treated as steady-state. If r>>ra does not hold, then the response of 
the measurement system is still relatively easy to analyze in terms of an average size 
extraction range over the sample period. If r is not large compared to r a> r x or r s then the 
problem becomes much more difficult. Particle diffusion makes theoretical analysis of the 
effect of transit time variation extremely difficult. 


DATA TRANSFER FUNCTION IN SCANNING MODE 

Consider a particle traversing the classifier from the aerosol entrance slit to the 
aerosol exit slit. It starts at radial position u= 1 on a streamline designated by Q in and 
ends at radial position w=7on a streamline designated by Q 0 ut- Q in and Q 0 »t are 
proportional to stream function. Q\ n is 0 on the center streamline of Q a , -1 on the 
downstream edge of the entrance slit and +1 on the upstream edge. Q ou t is 0 on the center 
streamline of Q s , -1 on the downstream edge of the exit slit and +1 on the upstream edge. 
In the steady-state mode, the dimensionless mobility Z p of a particle traversing from 
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(l,(jin) to (T.Qout) does not depend on the flow velocity profile. This is not the case in the 
scanning mode. 

Unless the flow profile can be accurately determined by numerical modeling, the 
DMA transfer function must be calibrated experimentally. The dimensionless flow velocity 
profile t/U is a function of position ( u,z ) and depends on the dimensionless parameters Re, 
(3, 8, 7, K a and Ks- Calibration is required over the entire range of these parameters 
anticipated to be encountered under ER -2 flight conditions. For a given DMA, 7, k, K a 
and Ks are fixed by geometry. /?, 6 are determined by the relative values of the DMA flows. 
If the flow split is kept constant then ( 3 , 8 remain constant. Re is likely to vary during 
flight due to the large pressure range covered. A constant mass flow control strategy would 
be required to keep Re constant. This is difficult in view of issues discussed below. 

If the functional time dependence of the scanning DMA voltage is of the form 
V— VVexp(f/r) where r is positive or negative, then all particles traversing the classifier 
from (l,Qin) to (7,Q 0 ut) follow the same trajectory. In this case, DMA transfer can be 
characterized by the function W= W(&n,Qout,S) where Sis the dimensionless scan rate. If 
Vis not of exponential form then the transfer characteristic depends on particle mobility 
(size) as well, making experimental calibration much more arduous. 

Assuming an exponential form for V, the sensitivity of the transfer characterization 
function W to the flow profile, hence, to Re, ( 3 , 8 , 7, k, /c a and /%, can be estimated by 
calculating IF for several trial profiles. The basic profile used consists of pure axial flow for 
the bulk of the region between u=l and u*=7 except for pure radial flow regions at the 
aerosol entrance and exit slits connecting streamlines in the entering and exiting aerosol 
flows to the corresponding streamlines in the axial flow region. Slug and Poiseuille axial 
profiles were used matched with uniform radial profiles in both cases. A uniform radial 
electric field was assumed. 

Particle trajectories were calculated for the above set of flow and electric fields. A 
particle entering on streamline at w=W2=l follows that streamline radially (at 
increased velocity due to radial electric field) to radial position u\n(Qm) where it enters the 
axial flow region. It moves axially due to the flow and radially due to the electric field 
across fluid streamlines until it reaches streamline Q 0 u t at radial position w 0 ut(< 2 out)- It 
then follows that streamline radially (at increased velocity) to radial position ^=^1=7 
where it exits the classification section. Wj n and Wout are determined by 
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Q( Win) = (1+Qin)^(l <5)/2(l+/?) 

QKut) = 1 - (l+Qout)/?(l+<5)/2(l+/?) 

where 

Q{u) = ■ f Mu') du ' 

is the fraction of total flow between radial position w and 1, and ^(w)=u z /C/is the 
dimensionless axial velocity profile. 


Wis the solution to the equation 

^in r. 


. 1 c (“Idu 

i = s-J u rfc) d ' 

The dimensionless mobility Z p of the particle traversing the classifier from (l.Qin) to 
( 7 j Q 0 ut) and located at radial position Wi n >w>w ou t in the axial flow region at time t is given 

by 


W = u + 


1-7 

LI+J3. 




If t 2> tin, t ou t and ti are the times the particle is at radial positions w 2 , Win, ^out and oj u 


e v 

e 

p (ir iout )/ r — — 


1 1 J / / — — 

- iV-Win 

- > (^Out~^in)/ 7 ’ _ 


1 

r h2 




^ 2 ” ^ in| 


/c a . 


L w J 


ir J 


1 

'hi] 


r/3(i+ffl 


t] - t out 

VV-Wout 

L «sj 


1 +P J 


Sr 


Derivations of these equations are given in the handwritten notes. 


The DMASCAN Fortran library is intended for use as an extension of the DMA 
Fortran library which is internally documented. DMASCAN contains the following 
functions with correspondence to the above notation as noted. 

FDENS(T,P) = p(T,p) - t It - ~ , x/ x 

ZSCAN(QSIN,QSOT,STIM,LPRO,IND,TSCAN) = [Z p e u ]{Qin,Qo\itA„{ti-t2)/T) 

EXPLIN(A,B) = root of e x +Bx-A 
WSCAN(WIN,WOT,GAM,STIM) = lF(wi n ,Wout,7,S) 

DSCANL(W.VF) = [%{u)-{W-u)/S\{uj,u z ) 

WFLOW(F,G) = u{Q) 

VFLOW(W,G) = • 
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The ZSCAN function returns the parameter Z v e i ^ r which can be used to calculate the 
mobility of a particle which traverses the indicated path and arrives at the exit slit (w=wi) 
at time t v LPRO is an integer flag indicating the axial flow profile: 0=slug, 2=Poiseuille. 
LPRO is passed to WFLOW and VFLOW through COMMON block FLOPRO. IND is an 
integer index referring to DMA configuration parameters stored in COMMON blocks 
defined by a call to DMASET in the DMA library. These COMMON blocks supply the 
required values of 7, *, P and 6 to ZSCAN. ZSCAN uses EXPLIN to solve the 

transcendental equations for the classifier transit times. The total dimensionless classifier 
transit time (i r i 2 )/r is returned as an output argument of ZSCAN. WSCAN finds W 
through a search routine which calls TSIMP, a Simpson numerical integration routine in 
the DMA library, for which DSCANL calculates the integrand. WFLOW supersedes the 
DMA library version for applications requiring LPRO=0. For the DMASCAN version to 
supersede the DMA version in linking, the DMASCAN library must be listed before the 
DMA library. FDENS is not used by any of the other DMASCAN functions but is used 
elsewhere for calculating Re. FDENS complements FVISC and FPATH in the DMA 
library which calculate air viscosity and mean free path. 


SCANF and SCAN are applications programs which utilize the DMASCAN library. 
Both programs are currently set up to give characteristics of the TSI Model 3071 DMA 
operated in the scanning mode at one atmosphere and 300K with 2 and 20 1pm aerosol ar.d 
sheath flows. Both routines ask the user for LPRO and 5. SCANF also requires and 

Q out and returns Wlni ofcut, V‘ /T and (irk)/ r - SCAN returnS ^ . and ( *„ r<a)/r 
f or Q 5n =o and Q 0 ut =0 and deviations of these parameters for an array of Q in and Qou t 

values specified by the user. The values are printed on the screen and into the file 
SCAN.DAT. Files of the form PnSx.DAT are sample results of SCAN output. 


The actual DMA transfer function ft as defined in Wang and Flagen (1990) is 
obtained by some sort of integration of Z v e^ T over Q in and Gout- This has not yet been 

derived. 

If mixing of the flow streamlines does not occur between the classifier aerosol exit 
slit and the CNC detection volume, then a rigorous analysis of the effective instantaneous 
DMA transfer function as seen at the CNC would require taking into account the varying 
transit time t d -U {U=U) between the DMA and the CNC for each exit streamline Q out . 
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p ARTICLE DIFFUSION 

Brownian diffusion of particles is an important effect for the proposed measurements 
of fine and ultrafine aerosols in the stratosphere. The particle diffusion coefficient is 
independent of pressure in the continuum regime and inversely proportional to pressure in 
the kinetic (free-molecule) regime. With a mean free path in air (at 23° C) of 0.067 fim at 
one atmosphere increasing to 1.36 fim at 50 mbar, the aerosol of interest is mostly in the 
kinetic and transition regimes. Consequently, particle diffusion is of increasing importance 
in the higher altitude range of the ER-2. 

There are two major effects of particle diffusion on the operation of the DMA on the 
ER-2. The first is particle loss to the walls of sampling and transport lines. For laminar 
flow these losses can be predicted using the theory of Gormley and Kennedy (1949). The 
second effect is a broadening of the response function of the DMA. This is due to diffusion 
within the actual classification section of the DMA and in the transport tube from the 
DMA aerosol exit slit to the detection volume of the CNC. This second effect is much 
more difficult to analyze and is also coupled with the loss effect. 

Cross-stream diffusion causes particles to travel randomly between adjacent 
streamlines in the classifier and transport tubes. For a particle entering a transport tube 
on a given streamline, this leads to a probability distribution of transit times and exit 
streamlines. The author is not aware of any analytic or even numeric or Monte Carlo 
solution to this problem. Particle loss to the walls of the transport tube should have some 
limiting effect on the breadth of this probability distribution. Streamwise diffusion 
contributes also contributes to the distribution of transit times. 


A rigorous analysis of diffusion effects in the classifier is even more difficult. 
However, Kousaka et al. (1985, 1986) and Stolzenburg (1988) have analyzed the problem of 
diffusion in the steady-state DMA. Kousaka et al. have done a dimensional analysis of the 
problem and show results of an approximate numerical model and experiments. 

Stolzenburg has constructed an approximate analytic model which is perhaps more 
convenient for estimating diffusion effects in the stratospheric DMA. 


The non-diffusing DMA transfer function H(Z P ) is triangular in shape (Knutson and 
Whitby, 1975a). As diffusion increases in importance, D becomes smaller and broader ar.d 
more Gaussian in shape. The analytic model of Stolzenburg shows that \JQ{Z p )dZ p ) 
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remains constant at «1 +S) while the mean (centroid) and standard deviation ot fi(Z p ) are 
given by 1+ff* 2 and [^(l+d 2 )/6+n- 2 (l+n- 2 )]i' 2 , respectively. The dimensionless d.ffusion 
parameter a- is given by c-'=G - £>. D- is the dimensionless diffusivity of the particle 
corresponding to Z p =l and is the same as the parameter D/S extracted in the dimensional 
analysis of Kousaka et al. (1985). G is a geometric flow factor evaluated as an integral 
along the particle trajectory. For the TSI Model 3071 with £=1/10, 6=0 and the same 
simplified flow Geld as described above, C=4.1104 for Poiseuille flow and C-3.5U0 for 

slug flow. 

The steady-state analysis of diffusion in the DMA can be used as a first 
approximation for the diffusion effects in the scanning DMA. As the relative scan rate S 
increases the error in this approximation also increases. The formulas above for the mean 
and standard deviation of Q{Z P ) can be used to estimate the degree of the diffusion effect. 
The full form of the steady-state transfer function for diffusing aerosols can obtained and 
employed in data reduction using the DMA library. 

The experimental data obtained by Stolzenburg (1988) to verify his diffusion theory 
indicate a discrepancy which may be significant for the stratospheric DMA. Figure 4.8 in 
Stolzenburg (1988) plots the relative error of the theory in predicting the mean mobility 
extracted by the second DMA in a tandem DMA (TDMA) experiment. Based on the 
experimental measurements, the theory shows an apparent growth of the aerosol in going 
from DMA1 to DMA2 when, in fact, there is no physical explanation for such growth. 

Other investigators have observed this same effect to a lesser degree for larger particle 
sizes. As seen from the top plot in figure 4.8, the effect increases as particle size decreases. 
The two experiments at different flow rates indicate that the effect also depends on 
something other than particle diameter. The lower plot shows a good correlation for both 
experiments with a measure of the relative broadening of the DMA transfer function due to 
diffusion. If this correlation is not merely coincidental it may indicate an error in the 
diffusion model which is proportional to the relative effect of diffusion. Depending on 
operational parameters, the relative diffusion effect in the stratospheric DMA is likely to 
extend beyond the range indicated in the lower plot of figure 4.8. Thus, substantial errors 
could occur in sizing the aerosol if the diffusion theory is in error as suggested. 

The experimental data of Stolzenburg indicate that the theory accurately predicts 
the width (standard deviation) of the transfer function. The discrepancy described above . 
may be due to an unpredicted skewing of the transfer function. There is some ambiguity in 
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the derivation of the form and use of G in the equations. The programs SCANG, GVAR 
and FVCOR represent attempts to estimate the magnitude of the effect of these 
ambiguities on the theory. They indicate that the magnitude is insufficient to explain the 
discrepancy in the TDMA experiment. The effect of space charge in the TDM A 
experiment was also analyzed and also appears to be insufficient to explain the discrepancy 
Hence, the discrepancy cannot be explained and could be significant for the stratospheric 

DMA. 


The DMASET subroutine in the DMA library evaluates G for the particle 
trajectory from Q in =0 to ~Q 0 ul =0 for Poiseuille flow. (Program GCON, not m the library, 
back calculates Gfrom variables in the DMA COMMON blocks.) The SCANG program, 
similar to SCAN, calculates G for an array of Q\ n and Q 0 ut values for either slug or 
Poiseuille flow. Contained within the SCAN.FOR file are functions FGEO and GINT. 
FGEO calculates G for specified Oi n and Q ou t values, flow profile flag (LPRO) and DMA 
configuration index (IND). FGEO calls GINT which supersedes the DMA library version 

for applications requiring LPRO=0. 



Appendix 2 


i 


AEROSOL CHARGING FOR STRATOSPHERIC 
nTFFFRHNTIAL MOBILITY ANAL\ ZER 


INTRODUCTION 

A major area of uncertainty in 

(DMA) operated m the stratosphere is the q aerosol in the sampled air, it is 

the aerosol. In order to obtain the size disUibut / n ), the fraction of 
necessary to know the form of the c arg ^ larity { This must be calculated for the 
particles of diameter D p carr} g g contr ibuted theories and/or experimental data 

GSmSSESSeasaiar™ 

theory works well in above. Hoppel and Marlow 

basis for most of the theoretical^ orKOi g p tQ obtain better agreement with 

have made the greatest efforts ^ 3 ^ he best data for making comparisons 
experimental data in the near ire pn ("\ qrO and of Pui, Fruin and 

appears to be that of Adachi Kousaka an uj ( conv enient manner. Adachi 
McMurry (1988) though the latter > 'S'^’^hTpui^pued his results to Fuchs’ 
compared his results to Fuchs' origma > ltheo y P^^ m PuPs paper , 

original theory and to a theon . j perreement with Pui’s data using a 

Marlow and McFarlane (1988) showed imp > _ ° e e k / 1989 ) showe d good agreement of 

modified theory. In another communication with Marlow) that 

their theory with Pui’s data and state ( p p ucbs > theory as modified by Hoppel 

there is an error in Marlow’s comment page . Thus ^ of th&t 

rYVara^e^eNn^nfL 1 ^ 

There are five major papers for 

foundation for calculating the lOMerasol attac t j p ] equations and shows some 
introduces (in his work) the t^e^ependent .ion J st comple te overall 

example solutions. Hoppel and Frick (198 ) g P reiterates most of the content of 
description of the aerosol charging p o ■ coefficients for typical tropospheric ions at 
Hoppel (1977) and includes tables solutions to the population 

normal temperature and pressur . aPrOSO i- a erosol interaction phenomena. Hoppel and 
balance equations including possible aero he ° a s t ° tachment coe ffi c ient theory against the 
Frick (1989) show probably the best te , r p uj p ru i n an d McMurry (1988). 

experimental ultrafine unipolar aer0 ?° 1 r c1 } 6 , J . d attac hment coefficients to ion 

Some information about the sen^tiu } lreat t ^ e time-dependent ion-aerosol 

properties is provided. Hoppel and Fr c ( ) ' h nonequ iiibrium nature of charge 

balance equations in depth, particularly different neutralizer 

distributions from neutralizers. Several solutions ; are shown ^™| e involved _ in 

configurations. The appendix proud T p„% r( \ to the uncertainty in 

calculating attachment coefficients particular y th . eg ^ thg three -body trapping 

determining ion properties. Note is made of an ff s t for -typical" neutralizer ions 

?a„.io„dd that there are many 
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typographical errors in the equations of these papers of which only a fraction have been 
discovered and noted. 


POPULATION balance e quations 


The Drst (las,) step in defining the aerosol charge ‘elates 

temS describing the birth and death 

processes for that species. 

death terms associated with each possible charge ate of an ^aeroso, p ^ ^ w 
also be death terms associated with losses to ' “ nta, "™ e neR atj v e'ions and free electrons), 
distinct ion species for a given polarity ’ is co ' 0 ’ di | additional birth/death terms 

StSSSSS* assumes there is only one 

distinct ion species of each polarity. 

There is one equation for each of the possible charge states ( 
particle The series o’f charge states is generally truncat^when ^P^obabdityrf 
achieving such a state becomes msigmCcan i . ^ ame size \ in adjacent charge states 

corresponding to ion attachment to part ( particles in the given charge state, 

and two death terms corresponding to , ion at modeled as a set of 

Aerosol coagulation is generally ignored. For ^a po yoisper e ; a ^ dass and the ion 

charge state and the ion-aerosol death terms in the ion 
equations become integrals over particular size. 

In addition to the time derivative equations then i may be &r 

Sed as totnl “vMch c^e the ion ime ^“^X^XlhS’Ste. 
respect to the ion-ion recombination/death term. Th rn n Cen trations In the case where 

Alternatively, the ion concentrations g , . zero'! then the aerosol equations 

ortli^rUo^oVVoP ve^ a^^bl'natiorr 

SMS S and tlJ ion equations need not be 

solved at all. 

In most ambient stratospheric s ‘‘ ua ^ eqi’atfons 0 

»*»•«>- Ifboth co " ditio " S h ° ld 3nd 1 
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neutralizer is not used ‘ e h d e ^^nij’ c ated above” TWs'maf not be the case for the High 

Speed Research P?o e gram. It depends o„ the rate and ^ *5SSta the 

exhaust. This also determmes « ma E “ ? des terns in the 
charfing T«Ss cannot t deUfned until the attachment coefficients have been 
determined. 

ATTACHMENT coefficients 

tonlero»UUachmet? coefficients^ required. The particle equations require only the .on- 
aerosol attachment coefficients. 

The ion birth term is a single number, q , indicating the ion generation rate per unit 
volume"' ^SK r «r.ti‘o« rate in the • ^ioactive 

and type, distance from the . S °J'“ roducfof the tonTrambination coefficient, a, and the 
recomb.nation/death term is J he prod recombination coefficient depends on 

CA»^fe^"!S2S5 

and mobility. The recombination rate must in the 
of generation rate and steady-state ion con balance equations need not be solved 

ion balance equations can be neglected, the f , V - S tate conditions hold, then 

but instead the ion concentrations may be mea and/or 

entire range of expected sampling conditions. 

The ion-aerosol birth/death terms are each a 

polarity i normalized to unit ion conce rharee level relative polarity and 

depends on ion mass and mob. ity and on . p H‘w® r 'f tocos 6 Besides indirect dependence 
dielectric constant (see e.j ., F ™» “ > “ d'reSiv fn temped 

^efficient is required to estimate the ion-aerosol trapp.ng d.stance. 

Thp t hprvrv for calculating an ion-aerosol attachment coefficient is rather 

v , j T , p Fnchs’ fl963”)°limitin^ sphere approach has been generally successfu . 1 
complicated. The huens ^lyooj nmiu 0 , P : cb hpre var j es from author to author 

si! 
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overall calculation arises. 

For pore coulombic attraction or repulsion the ion i» captured if 

trajectory takes it to the .actual parti c closest approach) for the possible ion 

impact parameter and minimum ap (P parameters approach the particle, pass 

trajectories. Ion trajectories l" h no apse 

through an apse and recede. Ion trajecto th^ i( P o{ &„ s captU red by the 

and ultimately spiral m to the particle su _ ■ numerical solution is required 

particle can be calculated analytically >" e.ther ca A f 

when there is both a net charge and an image char e^ ine ^ ^ herfc 

complicated if a distrib V qPJn Voppef (1977) assumed a Maxwellian velocity distribution. 
Following Keefe et al. (1?68), A j t j h e electric field of the particle. The resulting 

However, the ions are being hpt .,, ppn collisions with neutral gas molecules persists to 

velocity, bias built up by an ion between ^ 1W Krea ter mass of the ion. There may 

some extent after a single collision due to -the Jf 8^“ here which is not 

thus be an accumulated velocity bias of ^ J ons at the ^^nng^p ^ captured ion 
accounted for. Integration over p Ure cou l 0 mbic and image potentials with 

accuracy at reduced 

pressure and temperature is not easily extrapolated. 

In the limit of very small Particles the thl“he«y for S *on- 

approach the ion recombination coefficien . m , onarv and only the ion moves. In 

aerosol attachment generally treats the P ar 1 , t j ie j on ma ss and both must be 

the limit the particle mass approaches ; the same order in the th or , 

considered mobile. Hoppel (197 7) in Dr oblem is that the main mechanism for ion 

which tend to cancel out this error The other P r0 ^m s jnab opposite polarity 

Tecombination is three-body trapping. An ion approaching 

gains kinetic energy as it falls into the potential ^e this * xcess kinetic energy until its 

(presumably low probability), the ions rlptermine the fraction of ions which 

become ^ r ap p e d InVrbH^Jounfa'Pticle due to collision with one or more neutral air 
molecules. 

Hoppel references Natanson (I960) with 

electric forces limpact parameter less than the cntic^ . ^mpaa paramete The 

™ approximations concerning the velocity 

distribution and integration thereof. 

Of greater concern though is the definition of trapping distance. The excess energy 
from wSch it is defmedTs obtained from a similar delinit.on of the ron-ron trapp.ng 
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Since the analysis only looks at the of t otal energy which should be 

path length it appears that . ignores the of velocity 

distribution of ion impact parameters. 

4)lomnt . tn r .i c11 iate three-body trapping distances directly from theory have not 
Attempts to cal f ulat e Bru P e P kl f er (1964) as an example and uses his work 

the colliding ion and neutral molecule. 

on measurements of ion recombination coefficients. J 

Figure 3 in Hoppel (1977) and Figure 4 in Hoppel and Frich : (1986) ! nd !9J^ e ' t hat 

Snlnndi^aTnOTmll e^tion in Hoppe! and Frick 

HORfil I'd 12) gives dependence of trapping distance, a , on mean free path, A, 
rf 98 AlTl+5e2/3 6 fcrA)i'2 P -l] (attributed to Natanson). As pressure, p is reduced, goes up 
tomore sloily San A wU is inversely proportional to v The probab. htj -of collision 

- ite 

to the temperature de^ndences^uld be ofle^sermagmtu^de^han^ose^^t ^ three _ bo(Jy 

Snl i?mt be poSb"le lo ignore the effects of Lee-body trapping within the flight 
range of the ER-2 for particles larger than a few hundredths of a micron. 

himseXS 

kV l") d ' The°dS^ceme e „ C t U oFt HsTptolflomL'e ^(ide’LLe is'lessVn one mean 
f q - 16 A (1 % fA for A«a and s=A-a/3 for A»a where s is the displacement, A is tne 

mean free path' a/d a is the particle radius.) Hoppel 
attachment coefficient to this variation is unknown. 


TON PROPERTIES 

The hst steD in the aerosol charge distribution calculation is to obtain the necessary 
ion properties Th? properties needed are the mass, electrical mobi lily . d.ffusiv.ty mean 
ii-ormuf cnppd and mean free path. Using relations cited m Hoppel and Frick (1986) ano 
Pui, Fruin and McMurry (1988), the mean thermal speed can be calculated from mass an 



6 


the diffusivltv from mobility at a ;f. ve " 

ssissaa ^isSSi 

mass ratio indicates that this repre w p t n e D ersistence path might seem more 

A ^°A'SSS Snteffl to. the P Sbor KeXnard expressions by 
Sing the mass ratio factor. For m=200 am. and Af=29 amu the pers.stence path >s 
larger than the collision-free path by a factor of 2.8. 

Hoonel and Frick (1989) indicate that the calculation of attachment coefficient is 

SSSlSi 

various choices of the location of the Fuchs’ limiting sphere. A more careful sensiti y 
analysis is required to verify this. 

The minimum required ion properties are thus mass and . 1°' e xnJesSons tlVe 
and negative ions. Varions MpUta Hoppel 

and Frfck (IMoV in Adachi Kotsaka and Okayama (1985) and in Wiedensohler, 
LMkemeier, Feldpausch and Helsper (1986). It is not dear how moch confidence can be 
put in these expressions for the ER-2 application. 

Measurements indicate that a range of properties is to be expected for each ion 

TTnnnpl and Frick ('1990) Thus, ions produced in a neutralizer used on the ER - are 1 e y 
SSldSnS ambient stratospheric ions. In either case .t >s possibl ° 
?eat S d stinct Ion species individually in the population balance equa ions This would 

f s n ot recommended. This approach appears to be generally successful, but the error 
inoure^is^^ely to increase as the range of ion properties increases. and 

Fissan (1988) attempted to account for a second negative ion species - free elect ons - m 
their work with ultrapure gases. It seems unlikely that there will a ®iwufica 
concentration of free electrons in any charging situation imagined for an ER-2 based DMA 
but this should be considered. (See also O’Hara et al. , 1989.) 

The papers mentioned in this section provide a number of references which should 
nrovide a good stlrUng point for resolving the problem of determining ion properties. A 
few of these references S have been already obtained and are listed under ion references. 


CONnUlSlON’S 

Overall the problem of aerosol charging for a DMA operated in the stratosphere 
charging theory. 
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artifici^h^^^ 

concentration ratio in a bipolar g • , naturally occurring ones. The aerosol 

ssssffisSf '<% ^ ^ *““» o[ 

altitude than the naturally occurring charge distribution. 

If the naturally occurring charge <Ustr ^ t1 ^^ 
the only significant ion birth and death proc If t ^ ree -body ion trapping can be 

aerosol balance equations are relatively easy * b carried out with a reasonable 

required Zn^ZlTes “ dmoMite, The question remains as to what degree of 
confidence these may be predicted from the literature. 

The conditions stated in the previous pwagraph may “‘^fif^miltg 
measuring supersonic transport exhaust. T r > {j, e ion generation rate drop 

aerosols that have been m the d«l tor a long ^m et , o£th e 

distribution is likely to drop dramatically. 



Appendix 3 


Evaluation of Diffusion Batteries 
for Stratospheric Aerosol Measurements 


Diffusion batteries have the potential to be used to measure the size distributions of 
submicron particles in the stratosphere. This report presents the results of a feasibility 
study of this application of diffusion batteries. 

A diffusion battery is device used to separate particles on the basis of their diffusional 
mobility. Particles are pulled through a “battery” of screens or tubes. Particles of a given 
size have a known probability of striking the battery matrix. This probability is higher for 
smaller particles (with higher diffusion coefficients) than for larger particles (with lower 
diffusion coefficients). By varying the flowrate and/or the number of screens or tubes, 
particles of different sizes can be selected to penetrate through the battery to reach a 
counting device (usually a condensation nuclei counter). The controlling parameter in 
diffusion battery performance is the inverse Peclet number, <f>. For a tube-type diffusion 
battery, cj) is defined (e.g.. Brown et al., 1984) as 

<t> - N c 7tD*L c Q'\ 


where Nc is the number of tubes in the diffusion battery, D* the diffusion coefficient, Lc 
the length of the tubes, and Q the volumetric flowrate. Particle diameter enters the 
formula through the diffusion coefficient, given as 


kT 

37rnD P 


1 + 


2.512X 

D P 


0.87X. , 

— e xp ( 


-0.44D 

X 



where k is Boltzmann’s constant, T the temperature, r| the viscosity of the fluid (air), 1 the 
mean free path of the fluid and D p the diameter of the particle. The experimentally- 
determined efficiency of a particle penetrating a tube diffusion battery is given by (Brown 
et al., 1984) 

£=0.819exp(-3.650) + 0.097exp(-22.30) [<1>>0.04] and 

£=1-2. 560 273 + 1.20 + 0.177O 4/3 [0<0.04], 

Measurements are obtained from the diffusion battery by selecting a “stage”, a 
combination of flowrate, number of capillaries and capillary length. The number of 
particles penetrating each stage is then measured. These numbers (with error) and a 
knowledge of the penetration efficiency as a function of size (also with error due to the 
measurement of airflow) for each stage are used to extract the actual particle size 
distribution. The aerosol size distribution is assumed to remain constant while the number 



concentration of particles penetrating each stage of the battery is measured. If a single 
CNC is used to cycle through each stage, obtaining statistically reasonable counts may 
take several minutes to tens of minutes at stratospheric concentrations. 

A diffusion battery using commercially available glass capillary array filters has been 
modelled. The particle penetration efficiency for an ambient pressure of 50 mb at an 
instrument temperature of 273 K has been calculated as a function of particle diameter for 
each stage of the battery. These efficiencies are shown in Figure 2. The key feature of the 
penetration efficiency curves is the extremely broad nature of the response function for 
each stage of the battery. The measured response-the number concentration of particles 
penetrating each stage-is a convolution of the particle size distribution and these curves 
of particle penetration efficiency: 

Cj = J £j(D) N(D) dD, (3.1) 


where Cj is the concentration of particles penetrating stage j of the diffusion battery, the 
penetration efficiency of particles of size D penetrating stage j, N the concentration of 
particles of size D in the aerosol, and Di and D 2 the minimum and maximum particle 
diameters detectable with the CNC. Recovering the particle size distribution from the 
measurements Cj and known response functions Q(D) involves inverting the integral to 
recover N(D). Unfortunately, for any given set of measurements Cj, an infinity of possible 
solutions exist that satisfy(3.1). Additionally, because the problem is extremely ill-posed, 
small errors in measurements can lead to very large errors in the inverted size distribution 
using traditional matrix inversion techniques. For most applications, the greatest difficulty 
in using diffusion batteries involves this inversion process. Because this process is the 
limiting factor in using diffusion battery data for stratospheric aerosol size distribution 
measurements, considerable effort was expended in examining various techniques. 

The aerosol literature is full of techniques for solving the inversion problem. Most telling 
is that no single procedure has been universally accepted, and papers extolling new 
inversion techniques for diffusion batteries continue to appear. For this study, two new 
techniques for inverting diffusion battery data have been examined. The first of these, the 
Extreme Value Estimation (EVE) procedure, has been developed by P. Paatero of the 
University of Helsinki (Paatero, 1990). This technique does not attempt to find a single, 
“correct” size distribution, based on the argument that any size distribution that is 
consistent with the integral equation is valid. Instead, the EVE technique uses matrix 
inversion techniques to finds a solution set, limited to non-zero values, that is consistent 
with the measurements and their errors. The EVE program, which is commercially 
available from Dr. Paatero, was purchased for the purpose of this study. The source code 
for the software is not available. 

The EVE program was used to invert a number of modelled size distributions. For 
modelled stratospheric measurements, the errors in the measurements due solely to 
counting statistics were sufficiently large to render the results essentially meaningless. The 



deconvoluted family of size distributions that satisfied the equations within the estimated 
errors were broad enough to prevent reasonable determination of either the mean size or 
the geometric standard deviation of a simple lognormal, monomodal distribution. 

An alternative method, the maximum entropy (MAXENT) technique, is based on the 
argument that the “best” solution of the particle size distribution is that containing the 
least “information content” in the language of information theory (Yee, 1989). The 
solution with the least information content is the simplest and smoothest of the possible 
solutions; it has no extraneous information beyond that warranted by the measurements. 
This concept can be expressed mathematically by the Shannon-Jaynes entropy equation, 

H(n) = Jn(D)ln[n(D)]dD, 

Dl 

where H(n) is the Shannon-Jaynes entropy, and n(D) the normalized particle size 
distribution function (i.e., the particle probalility distribution) such thatJn(D)dD=l. 

Dl 

The MAXENT algorithm searches for n(D) that satisfies Eqn. 3.1 and that maximizes 
H(n). This method requires no a priori, ad hoc assumptions about the expected particle 
size distribution, as do most other inversion techniques, and produces a single size 
distribution by introducing the rational and valid constraint of maximum entropy. 

The MAXENT algorithm was obtained from Dr. E. Yee of the Canadian Defence 
Research Establishment. Again, the source code was not made available. As before, a 
number of modelled diffusion battery measurements based on a known particle size 
distribution function were used as input to the inversion algorithm. Measurement errors 
were also used as input, although the penetration functions were assumed to be perfectly 
known (no allowance was made for errors in flow measurement). The MAXENT 
algorithm showed excellent potential by correctly recovering some of the input size 
distributions with very good accuracy when measurement errors were small. However, 
the algorithm was inherently unstable; tests using realistic particle size distributions with 
realistic errors did not recover any of the input distributions and instead oscillated wildly. 
These failures may be in part a product of numeric difficulties in the complex algorithm to 
solve the maximum entropy equation. Until these limitations are reduced, the algorithm is 
of little use when realistic measurement errors are introduced. 

Conclusions 

Diffusion batteries coupled with CNCs show some potential for estimating the size 
distribution of stratospheric particles. The technique is severly limited by the lack of 
robust algorithms for retrieving data with measurement errors. The MAXENT method 
shows some potential when measurement errors are small. By artificially introducing small 
errors into the MAXENT program, some qualitative information on the mean particle size 
may be obtained, although the substantial errors in the details of the retrieved size 



distribution may result. Diffusion batteries cannot now be recommended for 
measurements of the stratospheric particle size distribution. 
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